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RECENT  DEVELOPMENT  AND  APPLICATIONS  OF  WENO 
SCHEMES 


CHI- WANG  SHU 

Division  of  Applied  Mathematics,  Brown  University 
Providence,  RI  02912,  USA.  E-mail:  shu@cfm.brown.edu 


Abstract.  This  paper  is  a summary  of  the  author’s  talk  given  in  the  Third 
AFOSR  International  Conference  on  DNS/LES  (TAICDL),  held  at  the  Uni- 
versity of  Texas  at  Arlington  on  August  5-9,  2001.  In  this  paper  we  briefly 
present  the  general  ideas  of  high  order  accurate  weighted  essentially  non- 
oscillatory  (WENO)  schemes,  and  describe  the  similarities  and  differences 
of  the  two  classes  of  WENO  schemes:  finite  volume  schemes  and  finite  differ- 
ence schemes.  We  also  briefly  mention  two  recent  developments  of  WENO 
schemes,  namely  a technique  to  treat  negative  linear  weights  and  a class  of 
high  order  central  WENO  schemes. 


1.  Introduction 

This  paper  is  a summary  of  my  talk  given  in  the  Third  AFOSR  International 
Conference  on  DNS/LES  (TAICDL),  held  at  the  University  of  Texas  at 
Arlington  on  August  5-9,  2001.  In  this  paper  we  briefly  present  the  general 
ideas  of  high  order  accurate  weighted  essentially  non-oscillatory  (WENO) 
schemes,  and  describe  the  similarities  and  differences  of  the  two  classes  of 
WENO  schemes:  finite  volume  schemes  and  finite  difference  schemes.  We 
also  briefly  mention  two  recent  developments  of  WENO  schemes,  namely  a 
technique  to  treat  negative  linear  weights  and  a class  of  high  order  central 
WENO  schemes. 

High  order  accurate  weighted  essentially  non-oscillatory  (WENO)  schemes 
have  been  developed  to  solve  a hyperbolic  conservation  law 

«t  + V-/(«)=0.  (1.1) 

The  first  WENO  scheme  was  constructed  in  [19]  for  a third  order  finite  vol- 
ume version  in  one  space  dimension.  In  [11],  third  and  fifth  order  finite  dif- 
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ference  WENO  schemes  in  multiple  space  dimensions  are  constructed,  with 
a general  framework  for  the  design  of  the  smoothness  indicators  and  non- 
linear weights.  Later,  second,  third  and  fourth  order  finite  volume  WENO 
schemes  for  2D  general  triangulation  have  been  developed  in  [5]  and  [9]. 
Very  high  order  finite  difference  WENO  schemes  (for  orders  between  7 and 
11)  have  been  developed  in  [lj.  Central  WENO  schemes  have  been  devel- 
oped in  [13],  [14],  [15]  and  [22].  A technique  to  treat  negative  linear  weights 
in  WENO  schemes  has  been  developed  in  [23].  In  this  conference,  Jiang, 
Shan  and  Liu  presented  their  new  results  on  developing  compact  WENO 
schemes. 

WENO  schemes  are  designed  based  on  the  successful  ENO  schemes 
in  [8,  26,  27].  Both  ENO  and  WENO  schemes  use  the  idea  of  adaptive 
stencils  in  the  reconstruction  procedure  based  on  the  local  smoothness  of 
the  numerical  solution  to  automatically  achieve  high  order  accuracy  and 
non-oscillatory  property  near  discontinuities.  ENO  uses  just  one  (optimal  in 
some  sense)  out  of  many  candidate  stencils  when  doing  the  reconstruction; 
while  WENO  uses  a convex  combination  of  all  the  candidate  stencils,  each 
being  assigned  a nonlinear  weight  which  depends  on  the  local  smoothness  of 
the  numerical  solution  based  on  that  stencil.  WENO  improves  upon  ENO 
in  robustness,  better  smoothness  of  fluxes,  better  steady  state  convergence, 
better  provable  convergence  properties,  and  more  efficiency.  For  a detailed 
review  of  ENO  and  WENO  schemes,  up  to  the  time  when  these  notes  were 
published,  we  refer  to  the  lecture  notes  [24,  25]. 

WENO  schemes  have  already  been  widely  used  in  applications.  Some  of 
the  examples  include  dynamical  response  of  a stellar  atmosphere  to  pressure 
perturbations  [4];  shock  vortex  interactions  and  other  gas  dynamics  prob- 
lems [6],  [7];  incompressible  flow  problems  [28];  Hamilton- Jacobi  equations 
[10];  magneto-hydrodynamics  [12];  underwater  blast-wave  focusing  [16];  the 
composite  schemes  and  shallow  water  equations  [17],  [18],  real  gas  compu- 
tations [20],  wave  propagation  using  Fey’s  method  of  transport  [21];  etc.  In 
this  conference  and  in  this  volume,  we  have  seen  several  new  developments 
and  applications  of  WENO  schemes  for  DNS/LES  of  turbulence  flows. 

The  organization  of  this  paper  is  as  follows.  In  section  2 we  discuss  two 
different  formulations  of  WENO  schemes,  namely  the  finite  volume  formu- 
lation and  the  finite  difference  formulation,  and  point  out  their  similarities 
and  differences,  both  in  one  dimension  and  in  multiple  space  dimensions.  In 
section  3 we  briefly  describe  two  recent  developments  of  WENO  schemes, 
namely  the  technique  to  treat  negative  linear  weights  [23]  and  the  very  high 
order  central  WENO  schemes  [22]. 
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2.  Finite  volume  and  finite  difference  WENO  schemes 

A finite  volume  scheme  for  a conservation  law  such  as  (1.1)  approximates 
an  integral  version  of  it.  Let’s  use  the  one  dimensional  example 

ut  + f(u)x  = 0 (2.1) 

to  illustrate  the  ideas.  Suppose  {/,  = , i]},  i = 1, N,  is  a parti- 

tion of  the  computational  domain,  and  A X{  = xi+\  — a^  i.  If  we  integrate 
the  PDE  (2.1)  in  the  cell  /*,  we  obtain 

[fHxi+i,t))  - /(«(®i_i,t))]  = 0 (2.2) 

where 

tt»(<)  = -r — / u{x,t)dx 

A Jit 

is  the  cell  average  of  u in  cell  7j . Notice  that  (2.2)  is  not  a scheme  yet,  rather 
it  is  still  an  exact  identity  satisfied  by  the  exact  solution  of  the  PDE  (2.1). 

This  is  the  starting  point,  though,  for  designing  a finite  volume  scheme. 

A semi-discrete  (discrete  in  the  spatial  variable  only)  finite  volume 
scheme  for  (2.1)  is  an  ODE  system  for  the  cell  averages  {u;(t)},  i = 1, ...,  N. 

In  order  to  obtain  such  a scheme,  we  would  need  the  following  reconstruc- 
tion procedure: 

Procedure  2.1:  Reconstruction.  Obtain  accurate  point  values  {u(xi+i,t)}, 
i = 0, ...,  N,  from  the  given  cell  averages  {u,;(f)},  i — 1, ...,  N. 

Notice  that  we  have  ignored  possible  problems  at  the  boundary.  We 
thus  assume  the  data  is  either  periodic  or  compactly  supported.  Boundary 
conditions  can  be  treated  in  a stable  and  high  order  fashion  depending  on 
the  type  of  boundary  conditions. 

WENO  is  simply  a specific  reconstruction  procedure.  Let  us  demon- 
strate the  fifth  order  version.  For  this  purpose,  the  approximation  of  {u(xi+ i , t)} 

uses  the  information  of  five  cell  averages,  from  the  stencil  {/,_  2,  li- 1,  /»,  li+i,  Ii+2)- 
This  stencil  is  not  symmetric  with  respect  to  the  point  sf+i  of  the  recon- 
struction. There  is  one  more  cell  to  the  left  than  to  the  right.  Thus  this 
reconstruction  is  good  for  upwinding.  The  procedure  consists  of  the  follow- 
ing steps: 

1.  We  break  the  final  stencil 

T—  {h-2,  li-i,  li,  li+i,  Ii+2}  (2-3) 

into  the  following  three  smaller  stencils: 

<5l  = {7i-2,7j_i,/j},  S2  = {Ii-l,Ii,Ii+l},  S3  = {Ii,  /j+l,  Ii+2}- 

Notice  that  each  small  stencil  contains  the  “target  cell”  Jj. 
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2.  We  construct  three  polynomials  pj  (x)  of  degree  at  most  two,  with  their 
cell  averages  agreeing  with  that  of  the  function  u in  the  three  cells  in 
each  stencil  Sj.  We  also  construct  a polynomial  Q(x)  of  degree  at  most 
four,  with  its  cell  averages  agreeing  with  that  of  the  function  u in  the 
five  cells  in  the  larger  stencil  T.  The  three  lower  order  approximations 
to  u(xi+i/2),  associated  with  Pj{x),  in  terms  of  the  given  cell  averages 
of  u,  are  given  by: 

1_  7_  11_ 

2 gUj_i  + g ttj, 

1_  5 _ 1_ 

~g^i-l  + gWl  + gWj-l-l,  (2.4) 

1_  5 _ 1_ 

-Ui  + -Uj+i  — QUi+2- 

The  coefficients  in  front  of  the  u could  be  derived  by  Lagrange  polyno- 
mials or  by  solving  a small  3x3  linear  system,  from  the  condition  that 
the  quadratic  polynomial  pj  (a;)  has  the  same  cell  averages  as  the  given 
u in  the  relevant  stencil.  See  [24]  for  details.  Each  of  the  Pj(xi+i/2) 
in  (2.4)  is  a third  order  approximation  to  u(xi+i/2)-  The  fifth  order 
approximation  to  u(xi+l/ 2),  associated  with  Q(x),  is  given  by: 

«fe+i/2)  = + 20»i+i  - 2o“i+2-  (2.5) 

3.  We  find  three  constants,  also  called  linear  weights, 

1 3 3 
71  =10'  72  =6’  73  = Io’  (2'6) 

such  that 


Pi  (*1+1/2)  = 
P2(*t+l/2)  = 

P3(*i+l/2)  = 


Q(%i+ 1/2)  = 71  Pl(*i+l/2)  + 72  P2(*t+l/2)  + 73  P3(*t+l/2) 

for  all  possible  given  data  iij,  j = i — 2,  i — 1,  *,  i + 1,*  + 2.  This  is 
to  say,  the  higher  order  reconstruction  Q(«t+ 1/2)  can  be  written  as 
a linear  combination  of  three  lower  order  reconstructions  Pj(®i+i/2)- 
The  linear  weights  given  in  (2.6)  depend  on  local  geometry  and  order 
of  accuracy,  but  not  on  Uj.  If  some  of  these  linear  weights  are  negative, 
special  techniques  must  be  used  and  will  be  described  in  next  section. 

4.  We  compute  the  smoothness  indicator,  denoted  by  f3j,  for  each  stencil 
Sj,  which  measures  how  smooth  the  function  Pj(x)  is  in  the  target  cell 
Ij.  The  smaller  this  smoothness  indicator  f3j,  the  smoother  the  function 
Pj( x)  is  in  the  target  cell.  In  most  of  the  currently  used  WENO  schemes 
the  following  smoothness  indicator  [11]  is  used: 


(2.7) 
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for  j = 1,2,3,  for  this  fifth  order  case.  These  smoothness  indicators 
are  quadratic  functions  of  the  cell  averages  in  the  stencil.  For  this  fifth 
order  case  they  can  be  worked  out  as  follows: 


A = 7“  (Ui- 2 - 2iq_i  + Ui )2  + \ (Ui- 2 - 4uj_i  + 3 Ui)2 


13 

12 

13 


A = yi  ( Ui -1  _ + wi+l)2  + \ (ui- 1 _ ^fi+1  )2 

13  i 

A = yy  (fij  _ + Ui+2 )2  + y (3fij  — 4fij+i  + 


(2.8) 


5.  We  compute  the  nonlinear  weights  based  on  the  smoothness  indicators: 


U)i 


UJj 


Zi*r 


UJj 


7; 


(e  + A) 


(2.9) 


where  7j  are  the  linear  weights  determined  in  step  3 above,  and  e is  a 
small  number  to  avoid  the  denominator  to  become  0.  Typically,  we  can 
take  e = 10~6  in  all  the  computations.  The  final  WENO  approximation 
or  reconstruction  is  then  given  by 


R{xi+ 1/2)  = wi  pi(a:i+1/2) +w2  P2(^+i/2) +w3  p3(^+1/2)  (2.10) 

With  this  WENO  reconstruction  procedure,  a finite  volume  WENO 
scheme  is  now  ready.  We  give  here  a very  brief  summary  of  all  the  steps 
of  a WENO  finite  volume  scheme  applied  to  (2.1)  in  the  one  dimensional 
scalar  case  for  the  positive  wind  f'{u ) > 0.  More  details  can  be  found  in, 
e.g.,  [11,  24,  25].  The  algorithm  consists  of  the  following  steps: 

1.  Given  the  cell  averages  ul  for  all  cells  for  time  level  n (starting  from 
time  level  0 which  is  the  initial  condition); 

2.  Reconstruct  the  point  values  A+i/2  ^or  ce^  boundaries  a?j+i/2  using 
the  reconstruction  procedure  detailed  above.  That  is,  we  use  (2.10) 
with  ujj  defined  by  (2.9),  using  7 j given  by  (2.6)  and  f3j  given  by  (2.8). 
The  superscript  ” in  j ,2  refers  to  the  fact  that  the  reconstruction 
has  a stencil  (2.3)  biased  to  the  left  relative  to  the  location  27+1/2-  This 
is  upwinding  according  to  the  wind  direction  }'{u)  > 0. 

3.  Form  the  residue  for  time  level  n in  the  method-of-lines  ODE 

~dt  ^ ~Kt,  (^+1/2)  “ /(“<- 1/2))  ’ 

and  move  to  time  level  n+1  by  a high  order  TVD  Runge-Kutta  method 
[26]. 
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Next  we  describe  the  setup  of  a finite  difference  scheme  for  solving  (2.1). 
A semi-discrete  finite  difference  scheme  for  (2.1)  is  an  ODE  system  for  the 
point  values  (uj(f)},  i = where  u;(f)  approximate  the  point  values 

of  the  solution  u(xi,t).  We  also  insist  on  a conservative  approximation  to 
the  derivative  f(u)x  in  the  form  of 

f(u)x\x=xi  ~ 2^1  (fi+l/2  ~ fi- 1/2)  (2-11) 

where  1/2  is  the  numerical  flux,  which  typically  is  a Lipschitz  continuous 
function  of  several  neighboring  values  of  Uj  around  X{. 

At  a first  glance,  the  finite  difference  scheme  has  nothing  in  common 
with  the  finite  volume  scheme  described  above,  as  they  approximate  differ- 
ent values  of  the  solution.  However,  the  following  observation,  first  intro- 
duced in  [27],  establishes  a close  relationship  between  the  two.  If  we  identify 
a function  h(x)  by 

/(«(*»  = nm,  (2.i2) 

where  we  have  suppressed  the  t dependency  of  the  function  as  we  are  in- 
terested now  only  at  spatial  discretizations,  then  by  just  taking  derivatives 
on  both  sides  of  the  above  equality  we  obtain 

/(»(*)),  = £ [*  (*  + if)  - h (x  - ¥).  • 

This  means  that  we  only  need  to  take  the  numerical  flux  as 

/i+1/2  = ^(^1+1/2)-  (2.13) 

If  we  could  get  an  approximation  to  h(x)  to  high  order  accuracy,  the  con- 
servative approximation  to  the  derivative  in  (2.11)  would  also  be  of  the 
same  high  order  of  accuracy.  Notice  that  (2.12)  can  be  written  as 

/(«*)  = hi, 

i.e.  we  are  given  the  cell  averages  of  h (since  we  know  the  point  values 
Ui,  hence  also  /(uj),  in  a finite  difference  scheme)  and  we  would  need  to 
reconstruct  its  point  values  h(xj+1/ 2)  for  the  numerical  flux  (2.13).  But  this 
is  exactly  the  same  reconstruction  problem  in  Procedure  2.1  above  for  finite 
volume  schemes! 

Thus  for  one  space  dimension,  the  finite  difference  and  finite  volume 
schemes  share  the  same  reconstruction  procedure,  applied  in  different  con- 
texts (on  cell  averages  of  u for  the  finite  volume  schemes,  and  on  point 
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values  of  f(u)  for  the  finite  difference  schemes).  They  involve  the  same 
complexity  and  cost.  Finite  difference  schemes  are  more  restrictive  in  the 
situations  that  they  can  be  applied,  as  they  only  work  for  uniform  or  smooth 
varying  meshes  and  a flux  splitting  (for  upwinding)  must  be  smooth. 

However,  for  multiple  space  dimensions,  there  are  essential  differences 
between  these  two  classes  of  methods,  when  the  order  of  accuracy  is  at  least 
three.  While  finite  difference  schemes  can  still  be  applied  in  a dimension- 
by-dimension  fashion  ( not  dimension  splitting!),  i.e.  computing  f(u)x  along 
a x-line  with  fixed  y using  the  procedure  above,  and  likewise  for  g(u)y , then 
adding  them  together  to  form  the  residue,  finite  volume  schemes  of  third 
order  or  higher  must  involve  multi-dimensional  reconstructions  from  cell 
averages  to  point  values  and  then  numerical  integrations  to  get  the  numeri- 
cal fluxes  along  the  boundaries  of  cells.  The  details  of  these  reconstructions 
can  be  found  in,  e.g.,  [24,  25,  9,  23].  As  such,  the  operation  count  and  CPU 
time  for  a finite  volume  scheme  is  around  two  to  four  times  more  expensive 
in  two  dimensions  and  around  five  to  nine  times  more  expensive  in  three 
dimensions,  compared  with  a finite  difference  scheme  of  the  same  order  of 
accuracy,  see,  e.g.  [3]  for  such  a comparison  for  ENO  schemes.  In  return,  the 
finite  volume  schemes  do  allow  more  flexibility  in  their  applications.  They 
can  be  applied  in  arbitrary  triangulations  and  do  not  require  smoothness 
of  the  meshes.  On  the  other  hand,  finite  difference  schemes  can  only  be 
applied  to  uniform  rectangular  or  smooth  curvilinear  coordinates. 

3.  Two  recent  developments  of  WENO  schemes 

In  this  section  we  briefly  describe  two  recent  developments  in  WENO 
schemes.  The  first  is  a technique  to  treat  negative  linear  weights  in  [23]. 
The  second  is  a class  of  high  order  central  WENO  schemes  in  [22]. 


3.1.  A TECHNIQUE  TO  TREAT  NEGATIVE  LINEAR  WEIGHTS 


As  we  can  see  from  the  previous  section,  a key  idea  in  the  WENO  recon- 
struction is  to  write  a high  order  reconstruction  as  a linear  combination 
of  several  lower  order  reconstructions.  The  combination  coefficients,  also 
called  linear  weights,  are  determined  by  local  geometry  and  order  of  ac- 
curacy. In  many  cases  these  linear  weights  are  positive,  such  as  in  (2.6). 
However  there  are  also  situations  where  some  of  them  become  negative. 
For  example,  if  we  have  exactly  the  same  setting  as  in  the  previous  section 
but  now  we  seek  the  reconstruction  not  at  the  cell  boundary  but  at  the  cell 
center  Xi , as  needed  by  the  central  schemes  with  staggered  grids  [13],  [22], 
then  step  1 would  stay  the  same  as  above;  step  2 would  produce 


1 _ 1 _ 23  _ 

— I/?'— 9 “f"  Uj  — 1 t~T  Uj  , 

24  12  24 


Pl(Xi) 
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(3.1) 


1 _ 13 _ 1 _ 

“ 24“i_1+12^  2AUi+U 

23  _ 1 _ 1 


24*  + 12^+1 


24 


Ui+2- 


Each  of  them  is  a third  order  reconstruction  to  u(xi).  The  fifth  order  re- 
construction to  u(xi),  associated  with  Q{x),  is  given  by: 


Q{xi)  = 640^-2 


29  _ 1067  29  _ 3 

Ji^-1  + 960  “ 480"i+1  + 640“i+2' 


Step  3 would  produce  the  following  weights: 


9 49  9 

80’  72  - 40’  73  “ 80’ 


(3.2) 


Notice  that  two  of  them  are  negative.  The  smoothness  indicators  in  step  4 
will  remain  the  same.  This  time,  the  WENO  approximation,  when  naively 
applied,  leads  to  unstable  results  because  of  the  negative  linear  weights.  As 
an  example,  in  Fig.  3.1,  left,  we  show  the  results  of  using  a fourth  order 
finite  volume  WENO  scheme  [9]  on  a non-uniform  triangular  mesh,  for 
solving  the  two  dimensional  Burgers  equation: 


in  the  domain  [—2,2]  x [—2,2]  with  an  initial  condition  uo(x,y)  = 0.3  + 
0.7 sin  (|(x  + y))  and  periodic  boundary  conditions.  We  can  see  that  serious 
oscillation  appears  near  the  shock  in  the  numerical  solution  once  the  shock 
has  developed.  The  oscillation  eventually  leads  to  instability  and  blowing 
up  of  the  numerical  solution  for  this  example. 

A simple  splitting  technique  of  treating  negative  weights  in  WENO 
schemes  is  developed  by  Shi,  Hu  and  Shu  in  [23]:  we  first  split  the  linear 
weights  into  two  groups 

7+  = \ (7*  + 3 |7i|)  , 7f  = it  ~ 7*7  * = 1,2, 3 

and  scale  them  by 


3= 1 


7 f = 7 f/a*, 


i = 1,2,3. 


For  the  simple  example  of  fifth  order  WENO  reconstruction  to  n(arj),  the 
split  linear  weights  are,  before  the  scaling, 

9 __  9 . 49  49  9 9 

7l  ~ 80’  7l  “ 40’  72  “ 20’  72  ~ 40’  73  _ 80’  73  “ 40' 
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Figure  3.1.  2D  Burgers’  equation.  Left:  fourth  order  WENO  result  at  t = 0.473, 
CFL=0.2,  without  any  special  treatment  for  the  negative  linear  weights;  Right:  fourth 
order  WENO  solution  at  t — 5/tt2,  CFL=0.2,  with  the  treatment  for  the  negative  linear 
weights. 


The  WENO  reconstruction  is  now  performed  on  each  group  separately,  by 
computing  the  nonlinear  weights  (2.9)  separately  for  with  the  same 
smoothness  indicators  /3j  in  (2.7).  The  final  WENO  reconstruction  is  then 
taken  as  a+  times  the  reconstruction  using  the  group  of  positive  weights 
minus  o~  times  the  reconstruction  using  the  group  of  negative  weights. 
The  key  idea  of  this  decomposition  is  to  make  sure  that  every  stencil  has 
a significant  representation  in  both  the  positive  and  the  negative  weight 
groups.  Within  each  group,  the  WENO  idea  of  redistributing  the  weights 
subject  to  a fixed  sum  according  to  the  smoothness  of  the  approximation 
is  still  followed  as  before.  For  more  details,  we  refer  to  [23]. 

We  notice  that,  as  the  most  expensive  part  of  the  WENO  procedure, 
namely  the  computation  of  the  smoothness  indicators  (2.7),  has  not  changed, 
the  extra  cost  of  this  positive/negative  weight  splitting  is  very  small.  How- 
ever this  simple  and  inexpensive  change  makes  a big  difference  to  the  com- 
putations. In  Fig.  3.1,  right,  we  show  the  result  of  the  approximation  to  the 
Burgers  equation,  now  using  WENO  schemes  with  this  splitting  treatment. 
We  can  see  clearly  that  the  results  are  now  as  good  as  one  would  get  from 
WENO  schemes  having  only  positive  linear  weights. 

3.2.  HIGH  ORDER  CENTRAL  WENO  SCHEMES 

In  [22]  a class  of  fourth  and  eighth  order  central  WENO  (CWENO)  schemes 
have  been  constructed  and  the  role  of  local  characteristic  decompositions 
on  eliminating  spurious  oscillations  is  demonstrated.  We  will  give  a brief 
description  of  such  schemes  here. 
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The  notations  are  the  same  as  before,  plus  the  new  one  for  staggered 
cells  Ii+ 1 = [xi,Xi+\}.  Let  At  be  the  time  step,  jxn=  u(xi,tn)  denotes 

the  point  values,  and  pn=  //,  u(x,  tn)dx,  1 = Si//  , u(x,tn)dx 

represent  the  cell  averages  at  time  tn  on  the  cells  It  and  /i+i,  respectively. 

The  CWENO  scheme  approximates  the  cell  averages  at  time  tn+1  based  on 
their  values  at  time  tn  with  staggered  space  grids.  We  integrate  (2.1)  over 
the  region  Ii+ 1 x [tn,  tn+1],  to  get  an  equivalent  formulation: 


vn+1 

*+5 


Ax 


rtn+1  „tn+ 1 

/ f{u(xi+i,t))dt-  f(u(xi,  t))dt 

Jtn  Jtn 


(3.4) 


What  we  want  to  do  is  to  find  approximations  of  the  cell  averages  un  i and 

1-\~  2 

the  two  integrals  in  (3.4).  Thus  the  algorithm  consists  of  two  major  steps 

to  evolve  from  {uf  to  { U71  x }: 

*+2 

1.  The  approximation  of  tt”  i from  the  knowledge  of  {uf  by  a WENO 

l+  2 

reconstruction.  Notice  that 


un 

Ui+ 


1 1 fxi+l  f%i+ 1 

J u(x,tn)dx  = J 7 u(x,tn)dx  + J u(x,tn)dx 


hence  we  only  need  to  reconstruct  ^ t u(x,tn)dx  for  all  i because 
— f ,+^  u(x,tn)dx  = uf  — f u(x,tn)dx  (3.5) 

l\X  J xi  /SX  Jx.  i 


by  conservation.  The  WENO  reconstruction  for  this  step  is  very  similar 
to  Procedure  2.1. 

2.  The  approximation  of  //"+  f(u(xi,  t))dt.  If  the  time  step  At  is  subject 
to  a restrictive  CFL  condition  At  < =f  - max  |//ti)|,  we  can  assume 
that  u(x{,  t)  is  smooth,  since  the  discontinuities  starting  at  tn  from  the 
staggered  grid  points  z,_i  and  x, ; , i have  not  reached  the  cell  boundary 
Xi  yet.  Hence  no  Riemann  solvers  are  needed  and  the  time  integrals 
can  be  evaluated  with  a quadrature  formula  to  high  order  accuracy. 
Notice  that  this  is  equivalent  to  a Lax-Friedrichs  scheme  and  the  same 
effect  can  also  be  achieved  without  a staggered  mesh  by  just  using  a 
Lax-Friedrichs  building  block,  such  as  those  WENO  finite  difference 
schemes  in  [11],  [1]  where  the  Lax-Friedrichs  building  blocks  are  used 
and  no  Riemann  solvers  are  needed  either.  We  could  for  example  use 
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a three  point  Gauss  quadrature,  obtaining 


f(u(xi,t))dt  « 


3 

At  ]T  aif{u(xi,  tn  + Air;)) 


i= l 


where  cq  and  tj  are  the  weights  and  knots  of  the  Gauss  quadrature, 
respectively.  Now  what  we  want  to  do  is  to  find  the  approximation  of 
the  point  values  u(xi,tn  + A/p)  from  the  cell  averages  {yf.  This  can 
be  achieved  by  solving  the  ODE  at  the  grid  points  x = xf 


u(xi,  tn)  « yn 


(3.6) 


with  |/7V:omputed  by  a WENO  reconstruction  from  the  cell  averages 
{ yf,  very  similar  to  Procedure  2.1  above.  The  ODE  (3.6)  is  solved  by 
a Runge-Kutta  method  to  obtain  the  approximation  of  u(xi,  tn  + Atri), 
with  the  aid  of  natural  continuous  extension  (NCE)  for  (3.6),  [29],  see 
[2]  and  [22]  for  details. 
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